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Thermal expansions and compressibilities are calculated from equations of state of the 
Griineisen form which are derived from a general partition function. The non-thermal energy 
and frequency are obtained as functions of volume, the constants in these functions being 
evaluated from data on the second virial coefficients of gases. Data of state are calculated for 
three types of thermal energy distributions, the classical, the Einstein, and the Debye. A 
limitation on the harmonic oscillator approximation is discussed as is the quantal effect in the 


law of corresponding states. 


HE harmonic oscillator model has long been 
used as a first approximation in theories of 
liquids and solids.! In this model monatomic 
molecules are assumed to vibrate independently 
with a Hooke’s law restoring force. The potential 
energy, U, consists of two terms, the potential 
energy of the molecules in their equilibrium posi- 
tions, Uo, and the potential energy due to the 
displacement, r, of the molecules from these 


positions. 
(1) 


where » and m are the average frequency of 
oscillation and the mass, respectively, U» and v 
being functions of the volume. 


U= Uo 


1 Among the more recent papers dealing with properties 
of state based on the harmonic oscillator model are: 
(a) K. F. Herzfeld and M. Goeppert-Mayer, Phys. Rev. 
46, 995 (1934); (b) L. Brillouin, Trans. Faraday Soc. 33, 
54 (1937) ; (c) J. F. Kincaid and H. Eyring, J. Phys. Chem. 
37 (1938); (d) M. Born, J. Chem. Phys. 7, 591 (1939); 
(e) G. Kane, J. Chem. Phys. 7, 603 (1939) ; (f) O. K. Rice, 
J. Am. Chem. Soc. 63, 3 (1941); (g) J. Walter and H. 
Eyring, J. Chem. Phys. 9, 393 (1941); (h) F. A. Matsen 
and G. M. Watson, J. Chem. Phys. 11, 343 (1943); (i) O. 
K. Rice, J. Chem. Phys. 12, 289 (1944); (j) M. M. Gow, 
Proc. Camb. Soc. 40, 151 (1944); (k) R. Fiirth, Proc. Roy. 
Soc. 183A, 87 (1944). 
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If the system is considered quantal the zero 
point energy, 3hv/2 must be included. Thus, U’ 
may be substituted for U defined above. 


(2) 


In this simple system the volume is a minimum 
at absolute zero, at which temperature the mole- 
cules are motionless if the system is classical, or 
possess the lowest amplitude of vibration if the 
system is quantal. As the temperature is in- 
creased the amplitude of the vibrations increases 
with the result that the volume increases. 

The Griineisen equation of state may be ob- 
tained from the following general partition func- 
tion for the solid state, 


Q=K(0/T) exp | — Uo/kT}. 


= 


(3) 


Here @ is the Debye temperature hy/k, and Up 
has the same significance as in Eq. (1). 


P=kT(6In Q/6V) r= —(6U0/6V) 
+kT(6In K(0/T)/5V)r. (4) 
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Fic. 1. Thermal expansion of neon, zero point 
energy neglected. 


Now 
kT(6 In K(0/T)/6V)7 
= —kT*(6 In K(6/7)/6T)y(6 In 6/6V)r 


= —E(T)(5 In 0/5V)r.? (S) 


On substitution of Eq. (5) the equation of state 
becomes 


P= 1n 0/6V)7; 


(8Uo/8V) E(T)(5 In ¥/V)r, 
since 


0=hv/k. (6) 


At zero pressure 


E(T)= In »/V) (7) 


d\n K(6/T) 


dln K(6/T)= d(6/T) (0/T). 


Then 


6T Vv d(0/T) \ 8T Jy 


_—dIn K(6/T) 8 
@(@/T) 
and 
T d(0/T) 6V Jr 
d\n K(0/T) 
8V /r. 


Then 


6T 


Inclusion of the zero point energy may be made 
by writing the partition function 


Q=K(0/T) Xexp [(—Uot+2hv)/kT], (8) 
which leads to | 
In 

—$h(av/aV)r. (9) 


At zero pressure 
(8U0/dV) r+ 3h(dv/dV)r 
(Ad In v/AV)r 
—(dU)/dV)r 
~ (an v/aV)r 


E(T)=-—- 


—hv. (10) 


This equation has been derived by several other 
methods.!*!f 4 

Three forms of the thermal energy, E(7), will 
be used :4 


VOLUME 


Z.PE. 
OMITTED 


TEMPERATURE - °K 


Fic. 2. Thermal expansion of neon, with and 
without zero point energy. 


3 E. Griineisen, Zustand des festen Korpers, Handbuch der 
Physik (Verlagsbuchhandlung, Julius Springer, Berlin, 
1926), Vol. X. See also Vol. XXIV/2. 

* Derived from the following partition functions: 

a. Classical: 
K(6/T) =(T/6)*; 


K(6/T) 


b. Einstein: 


c. Debye: 
T3 x3 
In K(6/T) =In (1—exp 
Op*J 
Op=hvm/k where vm is the maximum frequency in the 


Debye frequency spectrum. It is related to the average 
frequency, v, by the approximate relation v= jvm. 
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a. Classical 
E(T)=3kT; (11) 
b. Einstein 
E(T)=3kT 12 
(T)= (12) 
c. Debye 


6p/T x3 


E(T) =9kT(T/6p)! f (13) 


Both U, and » are functions of the volume; 
however, it is more convenient to use the distance 
between the equilibrium positions of the mole- 
cules, a, related to the volume by the following 


equation : 
V=Na?/y, 


where y is the packing factor equal to v2 for face- 


centered structures. 


CALCULATION 


In an earlier paper!® Uo and » were obtained as 
functions of a by averaging a Morse function for 
the potential energy of a pair of molecules, one 
molecule moving over the surface of a sphere of 
radius r, the center of the sphere being at the 
equilibrium position of that molecule while the 
other molecule is fixed a distance a away. When 
this average was multiplied by one-half the 
number of nearest neighbors, the quantity U was 
obtained. 


CA [—2n(a—R) ] 
4 


2na+1 
x| ( ) sinh 2nr—2n cosh 2nr| 
r 


—2 exp [—n(a—R,) ] 


na+1 
| ) sinh ur —n cosh nr|l. (14) 


r 


On expansion in terms of r this yields 
CA 
[—2n(a—R)] 
—2 exp [—n(a—R,) ]} 
CA 
exp [—2n(a—R) ](n’a—n) 


—exp [—n(a—R,) ](n’a—2n)}r?, (15) 
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Fic. 3. Thermal expansion of argon. 


the terms containing r* and higher being neg- 
lected. The use of the expanded form and the 
neglect of the higher order terms causes Eq. (15) 


to be a less exact description of the potential _ . 


energy for large values of a. This effect will be 
discussed later in connection with the calcula- 
tions on the thermal expansion and compressi- 
bility. Values of the constants A, m, and R, 
were obtained from second virial coefficient 
data via Lennard-Jones and Devonshire® and 
Buckingham.§ 

On comparison with Eq. (1) it may be seen that 


CA 
[2n(a—R) 


—2exp[n(a—R,)]}, (16) 
and 


exp [—2n(a—R,) ](n’a—n) 

—exp [—n(a—R,) ](n’a—2n)}. (17) 

Equation (17) may be expressed in terms of Uo. 

1d?U, 


1 


312 daz a da 


as the simple differentiation of Uo will show. 


5 J. E. Lennard-Jones and A. F, Devonshire, Proc. Roy. 
Soc. A163, 53 (1937); R. A. Buckingham, ibid. A168, 264 
(1938). Both of these references are quoted by R. H. 
Fowler and E. A. Guggenheim, Statistical Thermodynamics 
(The Macmillan Company, New York, 1939). 
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PRESSURE - ATM. 


VOLUME -CC 


Fic. 4. Dotted lines, classical behavior; full lines, 
quantal behavior. 


Originally the authors did not consider this to 
be a general equation but peculiar to the Morse 
form of U». However, it was later pointed out® 
that Eq. (18) can be derived for any analytic Up: 
See appendix. Since differentiation is simpler 
with the Lennard-Jones form than with the 
Morse, the former will be used. 


(19) 


where u and v are constants evaluated in refer- 
ence (5). Then 
dU, 3C 
= —(va*— 2), 
da 


v? = (5va*—22u), 


(20) 


(21) 
and 


da 2 da 
Finally from Eq. (7) 
6C (va* — 2u) (S5va*—22u) 
E(T)=— 
(40va* — 308u) 


1 (40va*— 308m) 


(Svas—22u) 


(22) 


qi? 


6 Private communication from K. F. Herzfeld. 
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and from Eq. (9) 
6C (va*— 2u)(5va*— 


E(T)=— 
(40va* — 308) 
3hp —C 
2 = 
The temperature corresponding to a given 
value of a may now be calculated for the various 
energy distributions. For the classical distribu- 
tion the right side of Eq. (23) divided by 3k gives 
the temperature directly. For the Einstein and 
Debye distributions the right side of Eqs. (22) 
and (24) divided by T and values of (E—Eo)/T 
obtained from the Landolt-Bornstein Tabellen for 
each distribution are plotted against 6/7.’ The 
Debye temperatures, 0p were obtained from the 
same tables, that for argon being 85° and for neon 
63°. The temperature was calculated from the 
value of 0/T at the intersection of the two curves. 


The results are plotted in Figs, 1 to 3. 
In Fig. 1 the thermal expansions resulting 


(Sva*— (24) 


PRESSURE -ATM. 


VOLUME -CC 


Fic. 5. Dotted lines, classical behavior; full lines, 
quantal behavior. 


7 @p is here regarded as an experimental constant chosen 
so as to fit experimental Cy data. Variation of @p is small 
compared to the variation of temperature so that 6/T 
may be regarded as a function of T only. See Fowler and 
Guggenheim, reference 5, pp. 146-149. 
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from the three energy distributions are compared. 
In Figs. (2) and (3) the zero point energy is 
included as formulated in Eq. (24). 

Equations (6) and (9), taking a and T as inde- 
pendent variables, yield Figs. 4 and 5. 


LIMITATIONS OF THE HARMONIC 
OSCILLATOR MODEL 


Rice,!' employing specific heat data to calcu- 
late a total energy from which he subtracts the 
thermal energy and zero point energy, has ob- 
tained the potential energy of a pair of argon 


molecules. This potential energy when plotted . 


against the distance between molecules has a 
very different appearance from the one which is 
obtained from second virial coefficient data. See 
Fig. 6. He states that this discrepancy may be 
caused by the deviation of the solid from true 
harmonic oscillator behavior, which is to be ex- 
pected only for small vibrations. He states further 
that even at absolute zero, the zero point vibra- 
tions may be large enough to cause serious devia- 
tion. In consequence it would seem that the 
foundation of the Debye theory and the Griineisen 
equation which are based on the harmonic 
oscillator model is less secure than is generally 
recognized. 

The results of this research are in accord with 
Rice’s view of the limited applicability of the 
harmonic oscillator model. Equation (14) for the 
potential energy of an oscillator takes the 
harmonic oscillator form only after expansion and 


cutting off the expansion at r?. Consequently this _ . 


formulation of the potential energy indicates that 
these systems can be harmonic oscillators only at 
small volumes and that the equations of state 
derived do not apply at larger volumes. This 
effect appears as an almost vertical thermal ex- 
pansion curve and a minimum in the compressi- 
bility curve at larger volumes. When zero point 
energy is included the effect is more pronounced 
as may be seen in Fig. 2 for neon, only a small 
portion of the curve being significant. For hydro- 
gen and helium the effect is even more pro- 
nounced, the zero point vibration being so great 
that the resultant large volume invalidates the 
harmonic oscillator formulation even at absolute 


zero. 
It is interesting to note that infinite thermal 
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expansion at larger volumes from the harmonic 
oscillator model have been obtained by Herzfeld 
and Goeppert-Mayer,'* Born,’ Kane,!* Gow,}! 
and Fiirth,!* employing more elaborate develop- 
ments. The first of these authors have considered 


° 2 4 6 \ 8 "A 
| 
NEON 
ARGON 
3 
@ 
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Fic. 6. Potential energy of a pair of molecules as a function 
of their intermolecular distance. 


the solid to be unstable at the point of infinite 
thermal expansion and have associated this point 
with the melting point of the solid, obtaining 
rather good agreement. If this were a legitimate 
association the harmonic oscillator could be re- 
garded a good model up to the melting point. 
Fiirth has examined this association in detail and 
has come to the conclusion that it is not justified, 
pointing out that solids would have high thermal 
expansions near the melting point on the basis of 
theory, when experimentally they do not.* The 
harmonic oscillator model, therefore, does not 
describe the behavior of solids at larger volumes. 


8 Fiirth points out as further argument against the asso- 
ciation that the ratio of calculated volume change from 
absolute zero to instability point to the experimental 
volume change from absolute zero to the melting point is 
always greater than one. The authors believe, however, 
that this is not a valid argument in the case of the inert 
gases since Fiirth, treating the problem classically, does 
not take into account the effect of zero point energy on 
the volume. See Fig. 2. 
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DISCUSSION OF RESULTS 


In Fig. (1) the thermal expansions resulting 
from the three energy distributions are compared. 
The Einstein form with 6=30p/4 is practically 
identical with the Debye form and need not be 
discussed separately. The classical form gives a 
higher thermal expansion than does the quantal. 

Classically at absolute zero one would expect 
the molecules to occupy positions corresponding 
to the minima in their interaction curves. See 
Fig. 6. If the molecules were vibrating at ab- 
solute zero by virtue of their zero point energy 
it is to be anticipated that the molecules would 
move farther apart as they do when thermal 
energy is added to the system. This effect is 
shown in Figs. 2 and 3. The effect is greater 
the lighter the molecule and as was pointed out 
above, is so large in the case of hydrogen and 
helium that the harmonic oscillator equation 
cannot be applied. This may explain in part why 
the experimental molar volumes of helium, hydro- 
gen, and neon are larger than one would expect 
from their interaction curves. 

In Figs. 4 and 5 the PVT curves are given, 
the classical solid being more compressible than 
the Debye solid. The deviation of quantal from 
classical behavior with regard to compressibility 
and thermal expansion is less for argon than for 
neon. 


APPLICATION TO THE LAW OF | 
CORRESPONDING STATES 

Neon, hydrogen and helium are abnormal 
liquids in that their properties deviate markedly 
from those predicted by the law of corresponding 
states, the deviation being commonly ascribed to 
quantal effects. Pitzer? has formulated the con- 
ditions which must be fulfilled to apply the law of 
corresponding states. The conditions relative to 
this discussion are that classical statistical me- 
chanics must apply and that the potential energy 
between pairs of molecules must be a universal 
function of R/R, where R is the distance between 
the molecules and R; is the position of minimum 
in the potential function. See Fig. 6. The 
statement is then made that properties of liquids 
should be compared at temperatures such that 
the ratios of the volume of the gas to the volume 
of the liquid are equal. When this is carried out 


°K. S, Pitzer, J, Chem. Phys. 7, 583 (1939), 
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neon, hydrogen and helium deviate markediy 
from normal behavior.!® 

For a given reduced temperature the ratio 
Veas/Viiquia Should be a constant for all normal 
liquids. However, for neon, hydrogen, and helium 
this volume ratio is lower than for normal liquids, 
indicating that the volume of the liquid is 
abnormally high. This is borne out qualitatively 
by the fact that in this research the quantal 
volume was found to be larger than the classical 
(i.e., normal behavior) at low temperatures. The 
effect was found to be smaller for argon than for 
neon. The difference between the classical and 
the quantal volumes becomes smaller at higher 
temperatures due to the greater thermal ex- 
pansion of the classical system, the quantal 
behavior becoming identical with the classical 
behavior at higher temperatures. Since argon is a 
liquid at higher temperatures and heavier than 
neon, these two effects combine to make argon a 
more normal liquid. By similar argument helium 
and hydrogen are less normal. 

The fact that quantal volumes are larger than 
classical volumes affects the comparison of other 
properties, though not markedly. Pitzer® has 
pointed out that this can account for only about 
1/7 of the total discrepancy (0.7 cal./mole/de- 
gree) in the entropy of vaporization of neon, the 
remaining difference coming from the quantal 
effect on the entropy itself." 


APPENDIX 
Let ¢(a) be the potential energy of a pair of molecules, 


3 
a distance a?= x;2 apart. The Taylor expansion ir. three 
i 


coordinates about a =a!, and neglecting cross differentials is 


From this expansion may be calculated the average 
potential energy ¢ of the pair of molecules for a displace- 
ment da from the reference distance a!. 


10 Quantal effects are important for light molecules and 
low comparison temperature, these two factors seeming to 
go together in that light molecules also have shallow 
interaction curves and consequently are liquids at low 
temperatures only. 

11'A very approximate calculation of the difference be- 
tween entropy of a quantal and that of a classical system 
may be made. At 20°K for neon, with @p=63 and 6=48, 


Spebye(9p) — Sciassical = Spebye(8p) — 3R—3R In (T/@) 
= 1.5 cal,/mole/degree 


= 


dp /Ax; = (0/Aa) (xi /a) 


and 


da? a? 
Replacing x;/a and x;?/a* by their average values 0, and 
4 respectively, 


_ +2(1-%5 do 


Then 


1d% 
2 da? 


a da 


dx;?, 


since the average of a sum is the sum of the individual 


¢=¢9(a') = 


SCATTERING OF LIGHT 


averages. Defining © 
dx? =da* as r*, 


iz) 


When ¢ is multiplied by } the number of nearest neighbors, 


C, there is obtained 
1d?Uo 


if2 1dU, 
3 dat ta |*. 
Now if U is expanded as a function of r 
U=U0+3(@U/dr*) or’, 
it will be seen that 


1/@U 


2 


i> 
~ 2L3\2 da?® 
which is identical with Eq. (18). 
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A molecular theory of the scattering of light by fluids of isotropic molecules is developed, 


utilizing recent advances in the statistical mechanics of condensed phases. The results are shown 
to be the same in first approximation as those of the continuous theory based on fluctuations 
originally proposed by Smoluchowski and Einstein. Higher approximations of the molecular 
theory are especially applicable to the study of critical opalescence and fluids consisting of very 


large molecules. 


INTRODUCTION 


HE theoretical calculation of the intensity 
of the light scattered by the molecules of a 
transparent substance is an old problem. Nearly 
half a century ago, Lord Rayleigh! developed the 
theory of the scattering of light by single 
isotropic molecules and extended it, with some 
success, to perfect gases. It was soon found, how- 
ever, that the intensity of scattering by con- 
densed phases was less than that predicted by the 
formula by more than one order of magnitude. 
This effect was rightly attributed to destructive 
interference between the wave trains scattered 
from different molecules, but the means of 
_ * Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 
** Present address, Brooklyn Polytechnic Institute, 


Brooklyn, New York. 
‘ Lord Rayleigh, Phil. Mag. 47, 375 (1899), 


calculating the extent of the interference were not 
available at that time. 

Subsequently the difficulty was elegantly cir- 
cumvented by Smoluchowski? and Einstein,* who: 
considered the liquid as a continuous medium 
troubled by small statistical fluctuations in 
density. The extent of these fluctuations could be 
calculated from the macroscopic compressibility 
of the medium, and the intensity of the scattered 
light was obtained without discussing the indi- 
vidual molecules at all. : 

It is now possible, however, asa result of recent 
advances in statistical mechanics, to obtain the 
Einstein-Smoluchowski formula from a molecular 
~ 2M. Smoluchowski, Ann. d. Phys. 25, 205 (1908). 

3 A. Einstein, Ann. d. Phys. 33, 1275 (1910). An error, 
which introduced a factor of (u?-+2)2/9 in the formula 


(24), is corrected by later authors. See references (7) 
and (8). 
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Fic. 1. Geometrical relations between light, system, 
and observation point. 


point of view. The treatment is also applicable in 
the neighborhood of the critical point, where the 
Einstein-Smoluchowski method encounters seri- 
ous difficulties. This molecular development is 
the subject of the present paper. 


FUNDAMENTAL EQUATIONS 


We consider a system of N identical molecules, 
which are smaller than the wave-length of light, 
in a compact volume V; N, V, and the tempera- 
ture, are such that the system consists wholly of 
one fluid phase. (See Fig. 1.) The system is ex- 
posed to a train of monochromatic, plane, and 
plane-polarized light waves whose electric field 
is given as a function of time and space by the 
equation : 


E(t, r) = Ep cos (ord), (1) 


where @ is a unit vector in the direction of 
propagation, A the wave-length in vacuum, c the 
velocity in vacuum, and ¢ the time. 

. In passing through the system, this wave train 
is modified to 


E’(t, r) cos (2) 


where uy is the refractive index of the system and 
E’, the “effective” field, is in general different 
from E in magnitude and direction. In each of the 
molecules of the system, the effective field induces 
a dipole moment p(é, r) : 


2 
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The polarizability tensor, a, is a molecular 
constant. We shall discuss in this paper only the 
case in which the molecules are isotropic ; that is, 
where @ may be adequately represented by a 
scalar. In addition, we shall not concern ourselves 
with variations in the value of Eo’ from point to 
point in the system, but shall assume it constant. 
The first restriction is for convenience only, to 
avoid an unduly cumbersome result ; the second 
is forced upon us by the present lack of a com- 
plete theory of the effective field E’. 

Under these restrictions p is the same as the 
average moment per molecule, p, of a part of the 
system containing many molecules, but in which 
the value of o-r does not change appreciably, 
which is known to be 


p(t, tr) =aE(t, r) 


cos —(@-ru—ct). (4) 
4nN 


Here &=[(u?—1) V/4aN ] is the average polariza- 
bility per molecule of the system to the external 
field, E. 

The light scattered by the system is observed 
at a point on the x-axis, x9, with the distance xo 
large compared to V?. Observation is made 
through an analyzer which only transmits light 
polarized in the z-direction. 

Consider a molecule, small compared to \, 
located within V at r;=ix;+jy:t+kz,. Ifi, j, andk 
are the usual unit ‘vectors in the x, y, and 2 
directions, the z-component of the electric field, 
measured at Xo, of the light scattered by this 
molecule is given by‘* 


m(Eo-k)(u2—1)V 
= cos —— 


x [u(o— i) ° ct |. (5) 


The z-component of the electric field, E,, of the 
total light scattered by all the molecules is the 
sum of their individual contributions: 


E,= (6) 
4 


4M. Born, Optik (Verlagsbuchhandlung, Julius Springer, 
Berlin, 933), Chap. VII, Sec. 81. 
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SCATTERING OF LIGHT 


The intensity of the scattered light, J,, is the 
time average (over an integral number of cycles) 
of the square of E,. The integral over time may 
be performed immediately; and the constant xo 
may be eliminated in the process: 


m(Eo-k)(u?—1)V771 
] 


| 


ost 


>’ cos 


i>j 


Qa 
a}. (7) 


The sums in (7) may be expressed and evaluated 
as integrals, provided certain distribution func- 
tions giving the probability of occurrence of 
specified values of the r; are known. The suc- 
cessful use of these functions is the kernel of this 
paper, and is possible only because of recent 
advances in the statistical mechanics of con- 
densed systems.* 

Mathematically, the argument is this: the sum 
of a function, f(r), over a large number of discrete 
values of the independent variable, ri, fo, -:-, 
etc. may be expressed as an integral over r 
provided the number of discrete values occurring 
in the interval between r and r+dr is known. 
This latter number is the probability of a value 
between r and r+dr, and is represented in the 
notation used by J. E. Mayer®® as (N/V) Fi(r). 


* A treatment of the problem in terms of a continuous 
rather than a molecular theory, which, however, makes 
use of the radial distribution function, has been given by 
— and F. Zernicke, Physik. Zeits. 27, 261 
a oi Mayer and E. Montroll, J. Chem. Phys. 9, 2 

$j. E. Mayer, J. Chem. Phys. 10, 629 (1942). 


A similar theorem holds for double sums: 


N(N-1) 


x f J f(t) f(ts)dridrs. (9) 


Here [N(N—1)/2V?]F2(ri, re) is the probability 
of the simultaneous occurrence of a value of both 
coordinates between rj, f2, and r;-+drj, 


DISTRIBUTION FUNCTIONS 


The distribution functions of the coor- 
dinates of the molecules, (N/V)Fi(r) and 
[N(N—1)/2V?]F.(r1, already defined in 
terms of probabilities, are so normalized that 
their average values in the coordinate space are 
respectively the number and a of pairs of 


molecules : 
J (-) F,(r)dr=N 
N(N—1) N(N-1) 
f J re)dridr, (11) 
2V? 2 


For fluids, where long range order is absent, 
F,(r) is a constant, equal to unity, and F2(r,, re) 
is a function only of the relative values of r; and 
fo, that is: 


(10) 


F,(t1, = F(t12), 


where is introduced. 

For the same reason, F2(r12) approaches a con- 
stant when fiz becomes large. Another useful 
function is then defined by 


= F.(f12) — lim 12). 


(12) 


(13) 


The integral, divided by twice the volume, of 
go(fi2) over the coordinates of the two molecules 
is a constant, written as be: 


1 


1 
f (14) 


The second half of the equation follows 
from (12). 
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EVALUATION OF I, 


In the present case, the f(r) of Eqs. (8) and (9) 
is sin (2xu/A)(e—i)-r. Using this, the sums in 
Eq. (7) may be changed to integrals ; 


1) 


N 
fre cos 


N(N-1) 
[fra T2) COS 


*T; COS (15) 


Since F;(r) =1, the first term is merely N/2, if 
edge effects be neglected. 


2m ON 


The evaluation of the second term is less 
trivial. Writing r; in terms of r; and rip so that 
we can write the 
integral as 


) J f F.(r12) cos 


Vv 


x (o—i) COS (17) 


This may be integrated over ri, if the second 
cosine is rewritten in terms of products of cosines 
and sines of r; and rie; the integral over r; of 
terms containing a cosine or sine of r; to the first 
power is zero, and to the second power is V/2, so 
that (17) becomes 


2. 
2V » 


A form more convenient for our purpose, 
is obtained by using (13) for F2(tiz). If, as usual 
we neglect edge effects,** the integral of 


** This neglect may be justified by. writing out the 
limits of both integrals in Eq. (16), and noting that the 
contribution to the integral of the surface is negligible. 


BRUNO H. ZIMM 


cos (2ru/d)(o—i)-fisdri2 over V is zero and the 
cross term finally becomes 


N(N-1) 


f cos —(e—i) (19) 
2V 

In the commonly met special case in which the 
Einstein-Smoluchowski treatment is useful, the 
integral (19) may be still further simplified. It 
may be seen from (13) that go(r12) becomes zero 
when ry, is large. Indeed, for ordinary fluids not 
in the vicinity of the critical point, it is well 
known that go(ris) is only different from zero at 
values of ri. less than several molecular diameters. 
If the molecules are small, cos (2ru/A)(#—i) -Ti2 
may be replaced by unity within this region, and 
(19) becomes simply b2/V: 


1 be 
(20) 


A value for 62, involving no significant ap- 
proximations, has been recently obtained by 
Mayer.® The derivation is too long to repeat here, 
but, from a generalized form of the grand parti- 
tion function, he was able to arrive directly at 


the value*** 
2V? 
2N2RT \ 


where z is the fugacity of the substance, nor- 
malized in terms of molecules per unit volume, 
and p is the pressure. By making use of the 
thermodynamic equation 


— (21) 


(0p/d2)r=NkT/zV, (22) 


we can transform (21) to the more usable form 


1 dln V V 
2 dp Jr N 


*** The F2(r12) used here differs from that employed by 
Mayer in reference (6), in that the latter defines his 
F.(f12) over an infinite system, whereas ours is defined 
over a finite system. In Mayer’s case, limris,.F2(f12) =1, 
but in our case for molecules of finite size, the limit differs 
from unity by a term of the order of magnitude of 1/N. 
Nevertheless, the go(ri2) and be as defined by (13) and 
(14) are identical with those used by Mayer. 


(23) 
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It should be noted that the only approximations 
involved in (21) and (23) are those inherent in the 
fundamental theorems of the semi-classical sta- 
tistical mechanics. ‘ 

If we now use (16) and (23) in (15), and 
overlook the truly negligible difference between 
N and N—1, we get the final expression for J, in 
the case of fluids of small molecules not near the 
critical point as 


I,= . 
2x07A4 \ ), (24) 


A useful related quantity is the ‘‘turbidity,”’ r, 
the extinction coefficient in Lambert’s law, 


T=Ige*'. (25) 


It may also be defined as the total intensity 
scattered per unit volume divided by the incident 
intensity, I> = Integration of (24) over all 
angles of the direction of the primary light (the 
details are given by Born‘) yields 


= 


This formula is in complete agreement with 
that derived by the usual Einstein-Smoluchowski 
methods.” § 


7S. Bhagavantam, The Scattering of Light and Raman 
Effect (Andhra University, Waltair, India, 1940). 

8J. Cabannes, La Diffusion Moléculaire de la Lumiére 
(Les Presses Universitaires de France, Paris, 1929). 


CONCLUSION 


In general, the approximation of replacing the 
cosine by unity in (19) will not be valid. In the 
neighborhood of the critical point, for example, 
or for fluids or glasses composed of very large 
molecules, go(f12) will not fall to zero at values of 
T12 less than the wave-length of light, and here 
use must be made of the complete form (19) 
rather than (20). The expression for J, then is 


2xA4 


V 
N 


The first term is just the light scattered by NV 
isolated molecules, while the second term de- 
scribes the interference resulting from correla- 
tions in position of the molecules. It will be noted 
that the second term is essentially the Fourier 
transform, as a function of the angles between 
illumination and observation, of go(ri2), in close 
analogy to the well known formulas for the 
scattering of X-rays. 

The extension of this method to many-com- 
ponent systems is straight-forward and will be 
reported in a later paper. 
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The paper discusses potentiometric measurements of hydrogen ion concentration carried out 
in hydrocarbon oils that have attained high electrical conductivity (a) through oxidation in the 
laboratory, (b) through deterioration in service, (c) through addition of certain compounds. The 
e.m.f. data are interpreted by comparing them with those obtained from lauryl sulfonic acid 
solutions in oil. The results indicate that the hydrogen ion concentration in the three types of 
oils is too low to account for the total observed electrical conductivity. It is obvious, then, that 
the predominant mechanism of ion generation is different from the usual acidic dissociation. A 
bimolecular reaction leading to two organic ions is suggested on the basis of evidence from 
organic chemistry and a quantitative relation that follows from conductivity data presented in 


this paper. 


1. THE PROBLEM 


N a previous publication! the author showed 
that hydrogen ion concentration can be meas- 
ured in dielectric liquids with a dielectric con- 
stant as low as 2.3, which is the usual order of 
magnitude in insulating liquids. The liquids 
previously investigated were xylene and dioxane. 
The author believes that, in the future, electro- 
.chemical measurements, not only with glass, but 
with other electrodes, will be as common in the 
realm of dielectric liquids as they are today in 
the domain of aqueous solutions. The purpose of 
the present paper is to show the application of 
the method to the field of electrical insulating 
oils, and to the investigation of their deteriora- 
tion in service. 

During oxidation of an insulating oil, its 
electrical conductivity increases, and this is one 
of the chief reasons that an oxidized oil is con- 
sidered objectionable. In the course of the 
oxidation various oxidation products, including 
acids, are formed. It is generally recognized that 
the increase in conductivity and the formation of 
acid are in some way related; in fact, the degree 
of oxidation is often estimated by the total 
acidity, or so-called acid number. While this 
point of view is certainly correct in a qualitative 
sense, the mechanism of this relationship has not 
yet been definitely established. 

One fact is certain: whenever the conductivity 
is independent of frequency (and this is very 
often the case), the conductivity is not due to 


1 Andrew Gemant, J. Chem. Phys. 12, 79 (1944). 


rotation of dipolar groups, but to the migration 
of ions. It is entirely unknown as yet, however, 
by what mechanism these ions are generated, and 
what their chemical constitution is. Referring to 
the statement made above, all one can say is 
that these ions are probably a product of the 
acids that are the result of oxidation. Lacking 
any positive knowledge in this matter, most 
writers working in this field assume, either im- 
plicitly or explicitly, that the acids formed (or, 
if metals were present during deterioration, their 
salts) are electrolytically dissociated. Among the 
authors stating this view, J. B. Whitehead,? J. C. 
Balsbaugh,* and the present writer‘ may be 
quoted. J. D. Piper® assumes electrolytic dis- 
sociation at the boundary of colloidal particle 
and solvent, and unequal distribution, between 
the solvent and particle, of the positive and 
negative ions formed. 

It was thought that use of the e.m.f. method, 
with glass electrodes, could throw some light on 
this problem. If dissociation of acids, in the 
usual sense, takes place in oxidized oils, there 
must be a hydrogen ion concentration corre- 
sponding to the electrical conductivity, since the 
latter would then be due to the presence of both 
organic anions and hydrogen ions. The method 
at hand allows a direct check as to whether or 


* J. B. Whitehead and F. E. Mauritz, Elec. Eng. 56, 
465 (1937). 

’J. C. Balsbaugh, A. A. Assaf, and J. L. Oncley, Ind. 
Eng. Chem. 34, 92 (1942). 

* Andrew Gemant, Liquid Dielectrics (John Wiley and 
Sons, Inc., New York, 1933), p. 61. 

Fam Sticher and J. D .Piper, Ind. Eng. Chem. 33, 
1567 (1941). 


146 


H-ION CONCENTRATION 


not this is the case. Whichever answer is ob- 
tained, it will help in narrowing the range in 
which the complete solution of the problem will 
be found. 


2. SOLUTIONS OF LAURYL SULFONIC ACID 
IN XYLENE 


In the solutions investigated during the course 
of the previous work on this subject, hydrocarbon 
oil was not used as the solvent. Before examining 
oxidized oils, it was necessary to study a system 
with an oil as the solvent in which a definite 
hydrogen ion concentration was established. 
Toluene sulfonic acid and picric acid used in the 
low molecular weight solvents are not soluble in 
oils. Piper* found that lauryl sulfonic acid is 
reasonably soluble in paraffin oil and gives conduc- 
tivities significantly higher than that of the oil. 
This acid was used in the present investigation. 
Because lauryl sulfonic acid as the solute was 
not included in the previous work, the first step 
was to examine solutions of it in a low molecular 
weight hydrocarbon, xylene, in order to estab- 
lish continuity with previous results. Because of 
the relative chemical instability of lauryl sulfonic 
acid, the results are less reproducible than those 
obtained with toluene sulfonic acid. 

Figure 1 shows the cell used for e.m.f. measure- 
ments which embodies several improvements 
over cells used previously in this study. A and B 


Fic. 1. Diagram of glass cell for e.m.f. measurements. 


are the two half-cells; A contains a standard 
solution of relatively high conductivity and 
density which fills the right-hand extension of A 


6J. D. Piper, A. G. Fleiger, C. C. Smith, and N. A. 
Kerstein, Ind. Eng. Chem. 34, 1505 (1942). 
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up to the opening F, B contains the solution to 
be measured. C is a large-bore stopcock, and D 
and E are the glass electrodes. The stopper that 
supports D is made of ground glass to fit the 
opening of A, the stopper for E is paraffinized 


DIELECTRIC 


LOGARITHM OF SPECIFIC CONDUCTIVITY, MHOS/CM 


-3.0 "25 -20 
LOGARITHM OF CONCENTRATION, MOLS/LITER 


Fic. 2. Specific conductivity vs. concentration for laury! 
sulfonic acid solutions. Solvent for Curves 1, 2, 3—xylene, 
for 4—hydrocarbon oil. 


cork. With stopcock C closed, A is completely 
(or nearly so) filled with the standard solution 
through stopcock G. The sample is introduced 
at the opening H; since the standard solution is 
made to have a higher density, no turbulent 
mixing will take place at F. During the measure- 
ment C is open; when the measurement is com- 
pleted, C is closed and the sample is drained out 
by using the three-way stopcock, J, and the cell 
is ready for another solution. After two or three 
measurements, the liquid surface of F is renewed 
by carefully opening C and G to allow a small 
amount of liquid to pass. 

For the standard solution an addition of about 
3 percent of ethyl alcohol and 6 percent of 
carbon tetrachloride to the solvent used in half- 
cell B is suitable. This increases the conductivity 
to the required level of 10-° mho/cm and also 
provides the necessary density differential. It is 
best to use the same solute as that present in B 
in a concentration of 0.01 or 0.02 normal. When 
the solvent of the sample was mineral oil, the 
necessary high conductivity could not be ob- 
tained; in this case xylene was used in the 
standard. Also, when oxidized oils were tested 
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DIELECTRIC 
CONSTANT 
1 2.75 
2 2.54 
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POTENTIAL, MV 
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Fic. 3. E.m.f. vs. conductivity ratio for lauryl sulfonic 
acid solutions in xylene. 


in B, the electrolytes in both half-cells were by 
necessity different; in this case toluene sulfonic 
acid or lauryl sulfonic acid was used in the 
standard half-cell. It is obvious that such a cell 
does not represent a concentration cell and would 
not be suited for a rigorous quantitative study. 
For the semi-quantitative problem at hand, the 
use of such a cell was deemed justified. 

The solutions investigated first were of lauryl 
sulfonic acid in xylene. By means of additions of 
ethyl alcohol, three solvents of dielectric con- 
stants 2.75, 2.54, and 2.38, variously, were ob- 
tained. The concentration of the acid ranged 
from 10-* to 10-* normal. The specific conduc- 
tivity is shown in Fig. 2 on a bilogarithmic scale, 
Curves 1, 2, and 3. The slope of the three curves 
is respectively 0.4, 0.7, and 1.3. From the mass 
action law it follows! that 


(1) 


in which = specific conductivity, Ao = equivalent 
conductivity, K =dissociation constant,c=molar 
concentration. This equation implies a slope of 
0.5. The gradual increase of the slope beyond 
the value 0.5 is in agreement with observations 


in the previous paper and also in papers of 
Fuoss and Kraus.’ 

The e.m.f. data are assembled in Fig. 3, the 
abscissa being the logarithm of the ratio of the 
conductivities in the two half-cells. For this 
case an approximate theoretical equation was 
shown! to hold, namely 


E= 118(1 —ty) logio (o4/on) 
+(1.55 (2) 


in which E=e.m.f. in mv, t7=hydrogen trans- 
ference number, o,4.and og=specific conductivi- 
ties of solutions in A and B, D4 and Dg=di- 
electric constants of solutions in A and B. The 
value of b is given by 1/b=1/r.-+1/ra, in which 
r,-=radius of cation, and rg=radius of anion in 
cm. The data for tz and 6 as computed from the 


experimental results are shown in Table I. The 


figures were obtained from the e.m.f. data; ty 
from the slopes and 6 from the intercepts of the 
curves in Fig. 3. The values for ty and b are 
both of the correct order of magnitude, showing 
a slight dependence on the dielectric constant. 

In these solutions, therefore, the e.m.f. is a 
function of conductivity, in exactly the same 
way as is to be expected from the reasonable 
assumption that the conductivity is due to a 
dissociation of the acid into hydrogen ion and 
anion. 

Before proceeding to test the same acid in 
hydrocarbon oil as solvent, it was necessary to 
decide in which way to take account of the 
rather different viscosities in the two half-cells. 
Because of the relation for the H+ ion concen- 
tration: 

[Ht ]=(10°/Ao)o, 


TABLE I. 


Ds tH b, cm 


2.75 0.46 * 
2.54 0.51 2.8X 10-8 
2.38 0.63 2.1X10-8 


* Since D4 was also 2.75, it was not possible to compute b in this case. 


the ratio of conductivities in Eq. (2) can be 
used only if the equivalent conductivities, Ao, 
are equal or nearly so. If this is not the case, the 


7R. M. Fuoss and C. Kraus, J. Am. Chem. Soc. 55, 
1019 (1933). 
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first term in Eq. (2) should read: 
118(1—tz) logio (4) 


or, if it is assumed that Ao is inversely propor- 
tional to the viscosity 7, and that the influence 
of other specific factors may be disregarded in 
first approximation : 


118(1—tz) logio (oana/onne). (5) 


Theoretically the term (5) should be prefer- 
able to the form given in Eq. (2), if the two 
viscosities differ. An attempt was made to decide 
experimentally between the two formulas by 
using the e.m.f. method. A solvent was prepared 
containing 6.7 percent polystyrene in xylene. 
This had-a viscosity at 25°C of 0.31 poise, 
whereas, that of the xylene was 0.0059. Two 
solutions of lauryl sulfonic acid (10-* and 3 X 10-*) 
were prepared in this solvent and the e.mf. 
measured against the standard solution of lauryl 
sulfonic acid. The data are plotted on Fig. 4 
which also contains Curve 3 of Fig. 3 for com- 
parison. The two points as plotted using the 
term (5) are indicated by @; plotted using 
Eq. (2), without viscosity correction, by O. 
The former points fit the extension of the line 
better than the latter. Thus, experimental evi- 
dence, in favor of the use of the viscosity cor- 
rection, was obtained. More work along similar 
lines should be carried out in order to clarify 
further this important question. This is, by the 
way, an example that shows the possible uses 
of the e.m.f. method in the field of liquid di- 
electrics. 


3. LAURYL SULFONIC ACID IN HYDROCARBON OIL 


In this section the results obtained on lauryl 
sulfonic acid solutions in a hydrocarbon oil will 
be reviewed. The mineral oil, which was a typical 
transformer oil, had the following characteristics 
(Table II). The oil was dried for a few hours at 
120°C before use. 


TABLE II. 


0.87 
0.14 
146 
— 50 
0.005 
1x10-" 
2.18 


Density, g/cm* 

Viscosity at 25°C, poise 

Flash point, C ' 

Pour point, C 

Acidity, mg KOH/g 

Elec. conductivity at 25°C, mho/cm 
Dielectric constant 
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The conductivity data are shown as Curve 4 
of Fig. 2: Its slope is 0.8, somewhat larger than 
that characteristic of typical dissociation, but its 
location with respect to the three curves referring 
to xylene solutions of higher dielectric constant 
shows that there is no break in the behavior of 
dielectric solutions down to a dielectric constant 
as low as 2.18. This is in keeping with the 
author’s idea of extending the e.m.f. method, so 
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Fic. 4. E.m.f. vs. conductivity ratio for lauryl sulfonic 
acid in xylene plus polystyrene. 


useful in highly conducting solutions, into the 
field of dielectrics. 

Data on e.m.f. are given in Fig. 5, the abscissae 
containing the viscosity correction. Dilutions 
from three different stock solutions are shown by 
different markings. The reproducibility of the 
data on these solutions is limited by the notice- 
able chemical instability of the acid, particularly 
in the presence of traces of moisture. In keeping 
with the general trend of this series of investiga- 
tions, no attempt was made to work under 
conditions of extreme dryness of the solutions 
tested. 

The intercept of the line with the ordinate 
axis is —45 mv. From this figure a value of 3, 
the average ionic radius, can be deduced, using 
Eq. (2) and taking ty=0.5 (the value for. the 
standard solution). D, is 2.6 and Dz is 2.2. 
Hence, }=12A units is obtained, five times as 
large as in xylene. Realizing that ions always 
carry molecules of solvation, the larger radius in 
oil, as compared with that in xylene, will be 
understandable. 
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Fic. 5. E.m.f. vs. conductivity ratio for lauryl sulfonic 
acid solutions in hydrocarbon oil. 


It should be emphasized, however, that the 
cell used for obtaining the data of Fig. 5 is not a 
concentration cell, since the solvents are quite 
different. The intercept, therefore, will contain a 
specific term, not covered by Eq. (2) which 
latter accounts only for the differences in di- 
electric constants. In spite of its limitations, the 
curve of Fig. 5 will serve as a basis for the 
interpretation of data on oxidized oils given in 
’ the next section. Figure 5 shows the e.m.f. vs. con- 
ductivity relation for a hydrocarbon oil in which 
' the total conductivity is essentially due to the 
presence of electrolytic dissociation products of 
an acid—hydrogen ions and corresponding anions. 


4. OXIDIZED HYDROCARBON OILS 


The results obtained up to this point permit 
an attack upon the problem outlined at the 
outset. In order to test oxidized mineral oils, 
five different samples were used. Four of these 
were transformer oils, deteriorated in service. 
The fifth sample was the same oil as that used 
in the preceding section, oxidized in the labora- 
tory at 150°C in the presence of copper as 
catalyst, and under free access of air. Table III 
gives the acidities of the five samples in question. 

Figure 6 shows the electrical conductivity as a 
function of concentration of the acid in mols/liter, 
as usual. The various dilutions were prepared by 
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Fic. 6. Specific conductivity vs. concentration of acid 
oxidized 


ydrocarbon oils. (Samples listed in Table 


means of the pure mineral oil mentioned before. 
The conductivities, measured with a.c. voltage 
ranging between frequencies of 100 and 10,000 
cycles, were independent of the frequency, thus 
indicating ionic conduction mechanism. 

Two conclusions can be drawn from Fig. 6: 

(1) The range of conductivities (10-" to 10-"), 
corresponding to molar acid concentrations be- 
tween 10-* and 10-, is about the same as that 
for the sulfonic acid solutions in oil (10-” to 
10-"). This appears to indicate a similarity 
between the two cases. 

(2) The slopes of the curves are respectively 
(starting from the top) 1.10, 1.05, 1.05, and 1.12, 
values that are close to 1.0. In contrast to this, 
the slopes for the strong acids were near 0.5 in 
the range of low concentration and decidedly 
above 1.0 in the range of higher concentrations, 
as the results in this and the previous publication 
show. Later in this paper this result will be 


TABLE III. 


Acidity in 
mg KOH/g 
0.42 
0.18 
0.30 
1.35 
0.57 


Number of 


oil sample Place of deterioration 


Transformer in service . 
Transformer in service 
Transformer in service 
Transformer in service 
Laboratory 
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tentatively interpreted from the standpoint of 
ion generation. 

The results of e.m.f. measurements on three of 
the five specimens are shown in Fig. 7. For 
comparison, the curve from Fig. 5 is also repro- 
duced. The following points should be men- 
tioned. 

(1) The potentials are not as well defined and 
stable as they were observed to be when sulfonic 
acid, or picric acid was present. They are un- 
stable and variable, of the type referred to as 
“floating” or undefined potentials. The unstable 
potentials are the more significant because the 
deteriorated oils were chemically rather stable, 
the conductivities not changing with time. Such 
floating potentials occur (also in aqueous solu- 
tion) when the concentration of ions correspond- 
ing to the electrode in use is so low that the 
actual ion concentration at the electrode surface 
depends on secondary uncontrollable factors. 

(2) The trend of the potential with varying 
dilution does not generally follow the direction 
characteristic of all previous measurements (in- 
crease of e.m.f. with increasing dilution). This 
again indicates that, if hydrogen ions are gener- 
ated at all by the acids present, their concentra- 
tion is too low to control the potential at the 
glass electrode. 
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Fic. 7. E.m.f. vs. conductivity ratio for oxidized hydro- 
carbon oils. (Solid line, referring to lauryl sulfonic acid 
solutions, reproduced from Fig. 5 for comparison.) 
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(3) The majority of measured points (eight 
out of ten) is considerably above the curve of 
comparison, which means that the hydrogen ion 
concentration is considerably below the total ion 
concentration. The amount of divergence can be 
estimated by reading from the extension of the 
comparison curve the abscissas corresponding to 
the measured potentials. It can be seen that 
conductivities due to hydrogen ions are about a 
hundred to ten thousand times smaller than the 
total conductivities. In other words, the hydro- 
gen ions—if present at all—are only a small 
fraction (10-* to 10-*) of the total ions present. 

These results are confirmed by a second set of 
data obtained with solvents of relatively low 
viscosity: a mixture of xylene and hydrocarbon 
oil. Two ratios were used, xylene-oil 1:1, and 
xylene-oil 2:3. Comparison was obtained with 
solutions of toluene sulfonic acid in these two 
solvents. The oxidation products were those con- 
tained in oil sample No. 1, the solutions being 
prepared by mixing oil No. 1, pure oil, and xylene 
in various required proportions. The results are 
shown in Fig. 8. 

These results are similar to those of the pre- 


XYLENE-OIL 


O TOLUENE SULFONIC ACID 
® OXIDIZED OIL No.! 


XYLENE-OIL 2:3 


OQ TOLUENE SULFONIC ACID 
@ OXIDIZED OIL No.! 


POTENTIAL, MV 


-0.5 0.5 10 
LOGARITHM OF CONDUCTIVITY RATIO 


Fic. 8. E.m.f. vs. conductivity ratio for mixtures of 
oxidized oil No. 1 and xylene. (Toluene sulfonic acid 
solutions in same solvent for comparison.) 
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Fic. 9. E.m.f. vs. conductivity ratio for solutions of 
phenol in xylene. (Solid lines refer to lauryl and toluene 
sulfonic acid solutions.) 


vious set. In particular it can be seen that the 
measured points are located well above the 
curve obtained by using toluene sulfonic acid as 
the electrolyte. In carrying out the semi-quanti- 
tative estimate as above, by shifting the points 
to the right until they reach the extension of the 
line below, again it will be seen that the hydrogen 
ion concentration is only a small fraction of the 
total ion content responsible for the high con- 
ductivity observed. 


5. SOLUTIONS OF PHENOL AND ASPHALT 


Two further sets of measurements, substanti- 
ating the main results obtained, are reported in 
this section. These two sets comprise e.m.f. data 
on solutions of phenol and asphalt in xylene and 
oil respectively. Phenol was selected because it 
is a typical weak acid in aqueous solutions 
(dissociation constant in water at 25°C, 1.3 
X10-) and, in this respect at least, represents 
the acids formed as oxidation products, which, 
as is known from titration in alcoholic solutions,® 
are also weak acids. Solutions of phenol in xylene 
have a conductivity appreciably higher than 
that of xylene alone. Asphalt was selected be- 


®R. N. Evans, and J. E, Davenport, Ind. Eng. Chem. 
(Anal. Ed.) 9, 321 (1937), 


cause it is known to cause high power factors in 
insulating oils.® 

The results of the tests involving phenol are 
shown in Fig. 9. One set (lower part of figure) 
refers to solutions in xylene (dielectric constant 
2.4) ; for comparison, the curve for lauryl sulfonic 
acid solution is given. The other set (upper part 
of figure) refers to xylene-ethyl alcohol (D=2.6) ; 
the comparison curve is that for solutions con- 
taining toluene sulfonic acid. The potentials ob- 
tained from phenolic solutions are not so irregular 
in their trend as those obtained from oxidized 
oils, but their position is above the lines of com- 
parison. If plotted upon the extension of the 
lines of comparison, they would indicate hydro- 
gen ion concentrations lower than that corre- 
sponding to their conductivity. 

The conductivity data for two different brands 
of asphalt in the oil characterized in Table II 
are given in Fig. 10, in which the concentrations 
are expressed in relative units. The slopes, par- 
ticularly toward lower concentrations, are close 
to 0.5. It is not clear, however, why these solu- 
tions should behave like typical weak electrolytes. 

The e.m.f. data for these two solutions are 
given in Fig. 11. Circles and crosses refer to 
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Fic. 10. Specific conductivity vs. concentration of two 
brands of asphalt in mineral oil. 


the same respective solutions as those in Fig. 10. 
The line of comparison is that for lauryl sulfonic 
acid. The potentials observed were of the floating 
type, and, with the exception of two points, are 
located above the line of comparison. The hydro- 
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H-ION CONCENTRATION 


gen ion concentration is certainly only a small 
fraction of the total ions present. 


6. DISCUSSION OF RESULTS. GENERATION OF 
IONS IN HYDROCARBONS 


The experimental results presented in this 
paper answer, in the negative sense, the question 
formulated in the introductory section. The 
hydrogen ion concentration, as measured by the 
e.m.f. method in oils possessing high conduc- 
tivities (as a result of deterioration in service, 
oxidation in the laboratory, or the addition of 
certain compounds), is too low to account for 
the total observed electrical conductivity. The 
conclusion is then inevitable, that the generation 
of ions from organic molecules present in oils of 
high conductivity does not take place by a 
process of typical acidic dissociation into H+ ions 
and an anion. If such a process is active at all, 
the acidic dissociation constant K of the com- 
pounds in question is too small to account for 
the major part of the ionic components present. 
In the previous publication! it was shown that 
the dissociation constant K of toluene sulfonic 
acid and picric acid in xylene and dioxane is of 
the order 10-?°—10-”. Hence, the K of the 
compounds in question must be several orders of 
magnitude smaller. The predominant type of 
dissociation in these solutions must be different 
from the acidic, and the majority of the positive 
ions different from H* ions. 

While a satisfactory answer to this latter 
problem can be found only by means of experi- 
ment, a guiding hypothesis will be developed 
briefly in the following. 

In aqueous and alcoholic solutions, acidic and 
basic dissociation is predominant, because both 
the H+ and OH™~ ions have a high energy of 
solvation. In hydrocarbons, however, acidic dis- 
sociation is exceptional, apparently being re- 
stricted to the strong organic sulfo and nitro 
acids. It is probable that, in such solvents, both 
the positive and the negative ion will be organic 
groups, and not simple ions such as H+ and OH-, 
because the affinity of an organic ion for the 
solvent molecules will evidently be greater, and 
the solvation energy sufficiently large, to com- 
pensate for the energy of dissociation. While 
organic anions are known in great number, 
organic cations are represented chiefly by the 
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various amines, and there is a great likelihood 
that it is this type of cations that is present in 
insulating liquids. Apart from the amines, it is 
known that organic cations exist in gaseous dis- 
charges ;* they were detected in large number by 
mass spectroscopy. As examples of such ions, 
C.H2+, CeHyt, CeHet, and the like, may be 
quoted. If such ions exist in a gas phase, their 
existence in a liquid dielectric is even more 
probable. 

The mechanism of ion formation of this type 
is not a dissociation, but probably a reaction of 
two or three molecules, accompanied by electron 
transfer from one to the other. While the usual 
dissociation is a monomolecular reaction, the type 
postulated here is a bimolecular or trimolecular 
one.* 

The conductivity data on oxidized oils (Fig. 6) 
supply experimental evidence in favor of this 
statement. If the reaction is A+B@Ct+D-, 
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then the mass action law gives 


(6) 


CaACB™Cc+Cp-, 


®G. Glockler, and S. C. Lind, The Electrochemistry of 
Gases (John Wiley and Sons, Inc., New York, 1939). 

* A continuation of this research now under way seems 
to substantiate these postulates. 
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c being respective concentrations. If the oxidized 
oil is present in a variable concentration, ¢o, in 
an indifferent oil, as was the case in the tests 
described above, then both ca and cg are pro- 
portional to co. As to cc+ and cp-, both are equal 
and proportional to o, the electrical conductivity. 
Hence, from (6) 
Co ~o", 
or 


(2) 


instead of Eq. (1), characteristic for dissociation. 
In logarithmic plotting, the slope in such a case 
must be unity, and, as may be recalled, this 
conclusion was well fulfilled in the data shown 
on Fig. 6. 

While the chemical nature of the ions C+ and 
D- could—as has been said—be found only by 
experiment, two known cases from organic chem- 
istry may be cited to illustrate that the formation 
of organic cations, apart from cases involving 
amines, is not purely hypothetical. 

One of the cases is that of cyclic hydrocarbons 
possessing double bonds, such as pinene and 
camphene. If an acid (usually a halogen acid) 
is added to CioHis, addition compounds like 
CioH 7X are formed. Meerwein has shown!® that 
these compounds have a tendency to form posi- 


como, 


10W. A. Waters, Physical Aspects of Organic Chemistry 
(D. Van Nostrand Company, Inc., New York, 1937). 
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tive ions CyHyt. Such ion formation may be 
assumed to take place to a small degree also in 
the presence of an organic acid HAc, thus: 


CyHis +HAc@CywH 


as an example of the above described general 
bimolecular reaction. 

As another illustration, triphenylcarbinol may 
be mentioned. As Gomberg and Anderson"! 
have shown, the latter is definitely basic in 
character, leading to the positive ion CPh,*. 
Such a dissociation has been known to take 
place in the presence of a halogen acid, but the 
possibility exists that it would also occur in the 
presence of organic acids; in this case the reaction 
would read: 


CPh;OH +HAc@CPhs3t +Ac~+H,0. 


While it is by no means maintained that the 
two types of compounds, as here discussed, 
actually play a role in the generation of ions in 
hydrocarbon oils, they are cited as evidence 
from organic chemistry in favor of the mechanism 
depicted above. 

The author wishes to thank his colleagues, 
A. G. Fleiger and N. A. Kerstein, for cooperation 
in this research. 


u1L. C. Anderson and W. A. Fisher, J. Am. Chem. Soc. 
66, 589 (1944). 
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The Calculation of the Contribution of Freely Rotating Groups to 
| Electron Scattering by Gases 
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The contribution from freely rotating atoms to the intensity of electron scattering by gaseous 
molecules is obtained by means of an asymptotic evaluation of an integral, using the saddle 
point method. It is possible now to estimate the effect of the approximation in using the Debye 
formula, which is a specialization of the derived series, and, if desired, to obtain a higher degree 


of accuracy. Two cases are considered numerically. 


N many molecules in the gaseous state internal 
vibrations of large amplitude and even prac- 
tically unhindered rotation are possible. This, of 
course, is well known from heat capacity and 
spectroscopic investigations. In his paper con- 
cerning the effect of various types of internal 
motion in gaseous molecules on the resulting 
intensity of scattered electrons, P. Debye! made 
an important contribution to the calculation of 
the effect of a freely rotating group in a special 
type of molecule by evaluating the integral 


sin sr 
{ —aw, (1) 
sr 


by means of an asymptotic formula. In the 
integrand, s=(42/X) sin (0/2), r is the varying 
distance of a rotating atom from another whose 
position in space is fixed, and which is not on the 
axis of rotation (for then r would not vary) and 
dW is the probability of finding the rotating atom 
in any position on the circle of rotation. s corre- 
sponds to sk in Debye’s treatment. 

It is the purpose of this paper to afford the 
investigator a means for estimating the effect of 


Fic. 1. 


1 P. Debye, J. Chem. Phys. 9, 55 (1941). 
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the approximation introduced by using the 
Debye formula and for obtaining a greater accu- 
racy, if desired. The general case, in which there 
are no restrictions on the molecular parameters, 
is treated and an expression is derived which is 
the sum of four asymptotic series. It can then be 
shown that Debye’s formula is obtained if one 
takes the first term in two of the series, neglecting 
the other two, and specializes it to Debye’s 
example. By means of these semi-convergent 
series it is easy to estimate the approximation 
introduced by using the Debye formula (or, more 
generally, the formula which is the sum of the 
first terms in two of the series) and, in addition, 
higher accuracy may be obtained by including 
additional terms. 


MATHEMATICAL DEVELOPMENT 


Consider the point, P, rotating with fixed 
radius of rotation, a, about the axis AC. Its 
distance, r, from a fixed point in space, Po, can be 
shown to be 


r=(c;—c2 cos B)', (2) 
where 


and is the angle of rotation. rmax=l2, fmin=h, 
and, clearly, 8=0 when r=]. 

The contribution to the total scattering func- 
tion from a pair of atoms at P and Po, respec- 
tively, is given by (sin sr)/sr. It is necessary 
however to consider the average contribution by 
the changing distance in the course of a rotation. 
Thus, we obtain an integral of the form of 
(1) which is an averaging of the contributions 
from all positions of orientation. We have by 
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expressing r as in (2) above, 


sin [x(1—p cos 8)! 3) 


x(1—p cos 


where x=sc;!, p=C2/c, and d8/2x=dW since, in 
free rotation, all positions are equally probable 
and the weight factor as a function of position 
reduces to unity. Integral (3) is the imaginary 
part of 


1 exp [ix(1—p cos 6)*] 


dp. 4 
ax Jo (1—p cos B)? 4) 


This integral will now be evaluated by the 
saddle point method, the details of which have 
been clearly presented by Debye.’ The principle 
of this method is that integrals of the type in 
(4) can be evaluated asymptotically by con- 
sidering the contribution in the vicinity of the 
extrema of f(8) = (1—p cos 8)!, which occur when 
B=0 and B=7, the saddle points. The evaluation 
is accomplished by choosing a contour of integra- 
tion along a path of steepest descent (steepest 
decrease of magnitude of the integrand) through 

the saddle points. 

In order to obtain now the contribution from 
the neighborhood of B=0, we define a new 
variable, 


t= f(B) — f(0) =(1—p cos (5) 


From (5) we get 
2(1—p)! (2—3p)t 
2p(1—p)! 


p 4p(1 


sin B(t) = 


2(1—p)! 
= 
pi n=0 


The values of bo, 61, be, and 63; are given along 
with the final expression. 

We now use (5) and (6) to express integral (4) 
as a function of ¢, and then make an additional 


2? P, Debye, Math. Ann. 67, 535 (1909). 
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change of variable, setting ixt= —7r. This gives 
exp [ix(1—p)*] 
ax(—ix)(1—p)ipi 
tb, be tbs 
x +++ (7) 


x x? x 


The choice of a contour of integration along the 
path of steepest descent through the saddle point 
allows us to choose ~ as the upper limit. Since 


f = T(v+1) 
0 


we obtain 


exp [ix(1— 


(1/2) 


be ibs 
4203/2) ‘| (8) 
x x? x8 


by choosing (—1)!=e7 ‘7/4, 

A similar term for the contribution from the 
neighborhood of B= can be obtained in the 
same way. By defining a new variable 


t= f(8) —f(x) =(1—p cos B)!—(1+p)4, (9) 
it is then found that 


2 i 
(1+?) (1 


(2+3p)ti 
sin B(t,) = 


2p(1+p)! 

ty? t;3 
(10) 
p 4p(i+p)} 

2(1 

p} n=0 
We finally obtain 
exp [ix(1+p)!—in/4] 


iby 
[ray 2)+—I'(3/2) 
x 


2 
x? x3 


FREE ROTATION 


We now have the complete evaluation of the 
original integral (3) in terms of the contributions 
from the two saddle points. It is the imaginary 
part of 


exp 
152b; 


ib, 3be 
x (1 + 
8x3 


Qe 4x? 
_exp 
3b.’ 
2x 


Ax? 8x3 


(12) 


and taking the imaginary part we have, 
1 
(wp) ix! 
X[S(x, p, + C(x, p, w/4)o2(p) 
+S(x, —p, 
+C(x, —p, —14/4)o2(—p) ], 
sin [x(1—p)!+7/4] 
(1—p)! 
cos [x(1—p)!+7/4] 
(1—p)* 


where 


S(x, p, m/4) 


C(x, p, /4) = 


2-3p _12—20p-+119? 
4p(1—p)} 32p?(1—p) 
407 1320 +166p*— 79" 
128p%(1—p)! 


The Debye formula is obtained if one expresses 
x and p in terms of s, J, and 1, sets o1=1 and 
o2=0 in (13), and if the molecule considered has 
the dimensions such that /,?—/,;?=4a?, where a is 
the radius of rotation. (An example is the case of 
a molecule with an identical group rotating at 
each end of a bond.) 


and 
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DISCUSSION 


It can be shown that the series in (13) are 
asymptotic in the sense of Poincare.* Thus these 
series are divergent and this property is mani- 
fested by an increase in magnitude of successive 
terms after a finite number of terms of steadily 
decreasing magnitude have been passed. In 
applying the series, the terms which are suc-. 
cessively increasing are not included and the 
approximation will be improved by taking more 
and more of the successively decreasing terms. 
The magnitude of the error resulting from cutting 
the series off at a particular term is roughly of the 
order of magnitude of the next term in the series. 

It is customary in using the visual method‘ for 
studying electron diffraction patterns of gas 
molecules to calculate theoretical intensity curves 
from s=2 to s=20 or further and to attach 
quantitative significance at s values greater than 
5. It is therefore desirable to evaluate the magni- 
tude of successive terms in the asymptotic series 
in the vicinity of s=5 for some typical molecules. 
(It should be noted however that with improve- 
ment in methods for obtaining accurate diffrac- 
tion intensities, it may be desirable to consider 
quantitatively the intensities at values of s less 
than 5.) 

In discussing the applicability of the derived 
series it is somewhat less cumbersome to refer to 
the complex expression (12). However, the re- 
marks to be made are readily applied to the final 
form (13). Two different molecules were chosen 
for discussion. Hexamethylethane was chosen 
because a study had been made on this molecule 
by S. H. Bauer and J. Y. Beach® using Debye’s 
formula and it represents a very favorable case. 
The second example concerns the rotation of the 
3 and 3’ carbon atoms in biphenyl, and, as will be 
seen, is less favorable for the application of the 
Debye formula. In Table I, the magnitude of the 
second and third terms included under the 
parentheses of each of the two series in (12) is 
recorded for three different values of s. The 
calculations were based on a hexamethylethane 


3E. T. Whittaker and G. N. Watson, Modern Analysis 
(The Macmillan Company, New York, 1943), American 
edition, p. 151. 

‘L. O. Brockway, Rev. Mod. Phys. 8, 231 (1936). 

5S. H. Bauer and J. Y. Beach, J. Am. Chem. Soc. 64, 
1142 (1942). 
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model in which the carbon to carbon distances 
were taken to be 1.54A and all the angles were 
considered tetrahedral. In Table II are recorded 
the same mathematical quantities for the 3 and 3’ 
carbon atoms in a model of biphenyl which has 


TABLE I. Evaluation of successive coefficients in the 
asymptotic series for a model of hexamethylethane. 


4x2 2x 4x2 


0.06 
0.01 
0.002 


| 3be2 by’ 3b/ 
2x 


0.13 
0.05 
0.03 


05 0.03 
0.005 
0.001 


s= 2 0. 
s= § 0.02 
s=10 0.01 


regular benzene rings with a carbon to carbon 
distance in the ring taken as 1.39A and a separa- 
tion between the rings of 1.52A. 

Whereas these two examples are by no means 
exhaustive, significant information may be ob- 
tained by examining the tables. The values for 
the terms included in the tables must be com- 
pared with unity, the leading term in the series, 
and thus the approximate magnitude of the error 
in omitting successive terms is readily obtained. 
For instance, for hexamethylethane the second 
term of the first series in expression (12) is only 
about 2 percent of the first term for s=5. The 
second term in the second series for s=2 is about 
13 percent of the first but this is in a region which 
is insignificant so far as the visual method is con- 
cerned. However for biphenyl in Table II, we see 
that the third term is larger than the second for 
s=2 and so calculations would include only the 
second term introducing a rather large error. It 
is also seen that rather large errors occur in 
omitting successive terms in the region of s=5 
to 10. 
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CONCLUSIONS 


Several conclusions can be drawn: 

1. Consistent with the spirit of the saddle point 
method, the accuracy of the asymptotic evalua- 
tion increases as s increases. 

2. The asymptotic series may be used to 
estimate the error involved in using the Debye 
formula and affords a scheme for reducing the 
error by utilizing additional terms. 

3. When the rotating atoms are far apart, the 
leading terms in the series may decrease in 
magnitude more slowly. This is illustrated in the 
case of biphenyl where the effect of the smaller 
value of » as compared to the value in hexa- 


TABLE II. Evaluation of successive coefficients in the 
asymptotic series for a model of biphenyl (re: 3,3’ carbon 
atoms). 


methylethane more than compensates for the 
effect of the increase in the value of x due to the 
larger separation. Since the barriers.which re- 
strict rotation are smaller when the rotating 
groups are farther apart, we can expect that the 
asymptotic series will find its greatest application 
in such cases. 
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25 4x2| 2x 4x2 
s= 2 — 0.23 0.28 0.28 0.32 
s= 5 0.09 0.05 0.11 0.05 
s=10 0.05 0.01 0.06 0.01 


